Underwater bubble pinch-off: transient stretching 
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At the point of pinch-oflF of an underwater air bubble, the speed of water rushing in diverges. Pre- 
vious studies that assumed radial flow throughout showed that the local axial shape is two smoothly 
connected, slender cones that transition very slowly (logarithmically) to a cylindrical segment. Our 
simulations show that even with initially radial flow, a transient vertical flow develops with com- 
parable speeds. Bernoulli pressure draws water into the singularity region while incompressibility 
forces it away from the neck minimum, generating significant vertical flows that rapidly slenderize 
and symmetrize the collapse region. This transition is due to a different mechanism, occurring much 
faster than previously expected. Vertical flows dictate the neck shape evolution. 



Introduction-yLathematical models of physical pro- 
cesses often predict the formation of a singularity. Ini- 
tially smoothly distributed and finite physical quantities, 
such as velocity and pressure, diverge in a finite time. 
Examples include supernovae, gravitational collapse into 
black holes, and the pinch-off of a fiuid drop. Most of 
the first-known singularities exhibited self-similarity and 
universality. Self-similarity indicates that the system, as 
it approaches the critical time (i — ^ t*) , is nearly iden- 
tical at different times except for some rescaling by a 
function of time to singularity {t* — t) , for example a 
system radius. Universality indicates that this function 
is independent of initial conditions and boundary con- 
ditions, so that every occurrence of a given type of sin- 
gularity occurs exactly the same way. An example of 
a universal singularity is the pinch-off of a water drop 
in air, in which surface tension smooths out azimuthal 
shape vibrations [H-Ill- We focus on the opposite phe- 
nomenon, namely the pinching off of an air bubble under- 
water (Fig.[T]). This commonplace phenomenon exhibits 
neither self-similarity nor universality, and its low energy 
makes it amenable to study using high-speed photogra- 
phy on a table-top setup p-Q- 

A nozzle submerged underwater blows a buoyant air 
bubble, which eventually detaches from the nozzle. Ini- 
tially, surface tension dominates and the bubble neck 
shape is quadratic, with a slight asymmetry due to the 
hydrostatic pressure gradient (Fig. [1^). As the water 
rushes in faster, inertial forces overtake surface tension 
as the dominant driver, while viscosity continues to be 
negligible. The neck shape takes the form of a hyperbola 
of rotation |9i]. As the system approaches the singular- 
ity, the two cones become more slender. Previous theory 
and simulation focused on the long-and-slender regime, 
where vertical flows are negligible a priori, so that indi- 
vidual vertical cross sections evolve independently (ol-fist . 
There, the approach to cylindrical is very slow (logarith- 
mic in i?rnin)- We focus on the transient regime where 
the vertical flows dominate the shape evolution and occur 
much faster. Finally, at small length scales, the singular- 



ity is pre-empted by one of two effects. One is airflow 
in the neck, which becomes important when the neck as- 
pect ratio reaches L^/Rmin ~ \/ puquid/pgas, or about 30 
for air- water systems [T5|. Second, if the neck is not 
perfectl y a xis y mmetric, azimuthal vibrations are gener- 
ated [7l.ll6l[T7|. As pinch-off proceeds, the amplitudes of 
the excited vibrations remain constant. When the neck 
radius shrinks down to the length scale of the vibrations, 
the sides contact before the void is filled. 
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Figure f : Pinch-off of bubble (dark region) from nozzle while 
submerged in water (white). Bright spots are optical artifacts, 
(a) Initially, the bubble neck shape near the minimum has a 
generic quadratic profile, (b-c) Near pinch-off, the shape be- 
comes two cones connected at the vertex by a short segment. 
A characteristic vertical length scale for the neck is the 
distance from the minimum to a height where the neck radius 
is 2_Rmin. (Images courtesy of N. C. Keim and S. R. Nagel). 

These vibrations are what excludes the bubble pinch- 
off from being self-similar and universal. Analogous in- 
stances of memory-encoding vibrations arise, for exam- 
ple, in the implosion of shock fronts jl^ and implosion of 
spherical voids. All such systems have two commonali- 
ties: damping forces become negligible compared to iner- 
tia, and the flow is predominantly inward collapse. The 
second condition, specifically radial fiow, was assumed in 
previous studies of azimuthal vibrations in bubble pinch- 
off fl^, Tt^I , so the behavior could be quite different with 
vertical flow present. 

In addition, we are interested in whether there is mem- 
ory of the axial shape. This type of memory is present 
when water pinches off in a bath of viscous oil: the axial 
curvature persists until pinch-off [l9|. In that system, the 
viscosity dampens any azimuthal vibrations. Since the 
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flow is purely radial, however, each cross section evolves 
independently. Incidentally, the collapse velocity is inde- 
pendent of height, and the curvature persists. 

Previous theory and simulation has only studied air 
bubble pinch-off in long-and-slender regimes, but pinch- 
off in other regimes is relevant. Inviscid collapse of ar- 
bitrarily shaped voids is important, for example, in cav- 
itation bubbles generated by complex flows, ultrasound, 
or laser pulses [20]; in hull breach; or in an object plung- 
ing into water [21|. The collapse of these voids is also 
controlled by inertia, but they are not necessarily in the 
long-and-slender limit, and significant vertical flows may 
be present. With that in mind, we simulate shapes with 
very squat and/or asymmetric cones under generic veloc- 
ity profiles. Remarkably, no matter how exaggerated the 
initial conditions, strong vertical fiows rapidly drive the 
shape to two slender, symmetric cones (assuming axisym- 
metry). This happens within observable time frames, be- 
fore the singularity is pre-empted. Therefore, every ax- 
isymmetric cavity pinches off with the well known long- 
and-slender behavior. 

Problem formulation-We use the fact that the inertia 
of the water inflow dominates viscosity, surface tension, 
and gravity in the limit the minimum radius i?min — ^ 0. 
Therefore, we set those ignorable terms to zero. The 
exterior velocity field u is incompressible (V • u = 0), ir- 
rotational (V x u = 0), and decays to zero far away. We 
also assume that the air in the bubble is dynamically pas- 
sive, with a uniform pressure Pair(^) whose value ensures 
constant bubble volume in time. 

Being curl-free, the velocity can be described by a 
scalar potential, defined as V0 = u. Equating the rel- 
evant stresses on the interface gives a non-linear differ- 
ential equation for 0, first-order in time and space: 



P. 



0. 



(1) 



The kinematic condition, which says that surface ele- 
ments are advected by ^ | ^ = V<p\g , determines how the 
interface evolves. Since the exterior flow is incompress- 
ible, V^c/) = 0, i.e., Laplace's equation holds. This allows 
u to be solved everywhere in the exterior given </> on the 
surface. 

We use the co-moving derivative — ^ + u ■ in 
equation [T] to evolve ip for discrete fluid elements. Also, 
we only need u± = Wcj) ■ n on the interface, so we use the 
Green's integral form for 4> on the surface: 



27r0 (xo) = 




(xi - Xo) 



xi - Xo 



|xi - xqI 



fids', (2) 




Figure 2: Pinch-ofT shifts and distorts the neck shape near 
the minimum from a severely asymmetric hyperbola into a 
symmetric one by reducing the opening angle of the larger 
"cone" until it equals the smaller opening angle, (a) Initial 
surface is comprised of a small bubble connected to a larger, 
elongated one. The upper opening angle 6+ — 40°, while the 
lower opening angle 6- =5°. The initial flow is radial, (b)- 
(c) successively magnify by 50x. Arrows indicate movement of 
the neck minimum, (b) At = 0.19, and (c) At = 4.4-10"*. All 
times non-dimensionalized by Ro/uq. Each innermost profile 
becomes the outermost in the next image. Every tenth node 
is shown. 



reduce to 1 dimension by assuming axisymmetry. At each 
time step, we use cubic splines to interpolate the N(t) 
discrete surface nodes x^. Then, for each i G {1 . . . N} 
we allow Xi to take the place of xo in equation [5] We 
use Gaussian quadrature to perform the integral over 
segments connecting adjacent spline midpoints, with the 
flrst and last segments being only half splines. This gives 
the needed relationship for {u.^} in terms of {ipi}. Us- 
ing 6t i?min/-Rmin, we advcct the N nodes using u and 
evolve (f) using equation [1] completing the cycle. 

In order to accurately resolve the pinch-off dynam- 
ics, we found it important to use a node distribution 
scheme that maintains a gradual variation in the spac- 
ing between node points and that continually adds nodes 
in the neighborhood of the minimum (see Fig. [5]). At 
each moment, if the vertical distance Az between nodes 
exceeds a maximum spacing Azinax(^j z) in the region 



where the neck radius is less than ARr, 



we add a 



new node at the midpoint of that spline. After ex- 
perimenting with several functions for Azmax, we found 
Azmax = max [0.15 min (|z| ,Lz) ,0.08Lz] allows the sim- 
ulation to accurately track the dynamics ^. 

We begin the simulation by specifying an interface 
shape and velocity held. In general, there are an infl- 



where Xq is a fixed point on S, Xi parameterizes the sur- 
face, and fi is the unit surface normal at xi pointing into 
the bubble 

With both governing equations defined on the inter- 
face, the formulation reduces to 2 dimensions. We further 



^ For the up-down asymmetric case, we define the current height of 



the neck minimum as zq, R \^zo ± j = 2/?niini and Azniax = 

max 1^0.15 min — zo| I ^z"^ I 0.08 min J , L j for z > zq 
(reverse superscripts for z < zq). 



Figure 3: Pinch-off stretches a bubble neck that initially is symmetric about the minimum toward a cylindrical segment, (a) 
Hyperbola with opening angle 9 — 75° and radial flow, (b)-(d) successively magnify by 12x. Time between snapshots: (b) 
At = 0.28, (c) At = 2.4 ■ 10"'', and (d) At — 7.4 • 10~®. Each innermost profile becomes the outermost in the next image. 



nite number of each. Here, we exaraine generic smooth 
distributions. To expedite the computation, we do not 
prescribe a purely quadratic shape profile at t — 0, but 
instead use the result from previous studies that a slender 
quadratic neck evolves into a hyperbola, and prescribe: 



R{z) = Rc 



\ 



1 + tan-' e I — 
Ho 



(tan2 e) {z^y 



A{Rl 



where 29 is the opening angle of the cone, Rq = 
Rmin = 0) , and i?o = 10i?o is the radius of the 
end cap. We experimented with different smoothly- 
varying initial velocity fields, with corrections in the 
end caps to preserve bubble volume. One such field 
is given by the normal velocity on the surface: uj_ = 



7S*. fncompress- 



Mo^o [i?o + tan2 e (l + tan2 ^) z2] 
ibility then determines the tangential component of the 
velocity. The 75^* is a correction term necessary for bub- 
ble volume conservation, where s is the arc length along 
the surface from the neck minimum and 7 is a con- 
stant. The parameter ^ specifies the initial orientation 
of the vertical velocity in the neck region. If ^ > 6*, 
the vertical velocities are directed toward the midline, 
compressing the neck aspect ratio into a more squat 
shape. This configuration, though, is unstable. Very 
soon (i?min ^ 0.1i?o) the vertical velocities near the mid- 
line flip to orient away from the midline and proceed to 
stretch the neck. If, on the other hand, the system is 
initialized with ^ < 0, the vertical velocities point away 
from the midline and continue to point away, stretching 
the neck for the entire collapse. After their respective 
transients, however, both situations fall into the same 
dynamics, only differing by a time offset. Therefore, for 
the remainder of the paper we consider an initial radial 
flow, which corresponds to ^ = 6*. To specify a radial flow 
for up-down asymmetric shapes, we use an explicit radial 
flow with end-cap corrections: uj_ = ~""'^° f • fi — 7s'*. 

Results-Small values of 9 with radial flow correspond 
to the regime of the slender-body approximation em- 
ployed in previous works. Our results (not shown here) 
agree quantitatively with their results in this regime. 
Here, we focus on large values of 6. Since the bubble 
break-up experiment shows a slight up-down asymmetry, 
we first examine a severely up-down asymmetric shape. 
We produce such an initial state (Fig. 2a) by using differ- 




Figure 4: (Color online) A toroidal pressure peak develops 
just inside the neck minimum that has a cross-section aspect 
ratio roughly Lz by Rmin and tracks the neck inward with 
speed i?min. Water is accelerated away from the peak, gen- 
erating vertical velocity on the surface comparable in magni- 
tude to radial velocity. The point of max on the surface 
is indicated. This pressure distribution is robust to changes 
in initial shape and velocity and develops quickly. The case 
shown is Rmin/Ro = 0.034 for an initial shape with 9 — 75° 
and radial flow. (Top) corresponds to Fig. [3](b), innermost 
profile. (Bottom) zoomed 5x. Pressure is in units pug. 



ent values of 9 for the top and bottom portions of the bub- 
ble instead of prescribing a hydrostatic pressure gradient. 
The upper portion of the bubble is chosen to have the 
larger opening angle. We set the initial flow to be radial 
in the neck region. Figs. [5Jd-c show the shape evolution. 
Initially, the minimum shifts toward the smaller cone. 
As the minimum radius decreases, the upper portion of 
the bubble surface rolls while the lower portion more or 
less preserves the same profile. As a consequence, the 
profile near the minimum quickly approaches a symmet- 
ric shape. Simulations starting with various choices for 
up-down asymmetry yield the same qualitative outcome 
(not shown). We have found that an up-down asym- 
metric initial state always proceeds through a 3-stage 
evolution. First, it evolves into a symmetric profile as 
described. Second, the symmetric cones become slender 
near the pinch-off point. Finally, there is a slow transi- 
tion to a cylinder. Hereafter we will focus on the second 
stage, which connects the asymmetry in experiment to 
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prior works focusing on the last stage 
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Figure 5: (Color online). (Inset) Both the maximum vertical 
velocity on the surface, u^, and the radial collapse speed -Rmin 
diverge as pinch-oft approaches. (Main figure) Dividing the 
vertical velocity profiles from different times by 7?min and z 
by Lz scale out most of the variation, showing that the veloc- 
ity field is controlled by the break-up dynamics. The rescaled 
profiles evolve non-monotonically. The vertical velocity ra- 
tio at early times (I, -Rmin = 0.9-Ro) reaches a maximum at 
-Rmin ~ 0.075-Ro, (II), then broadens in space (III). Eventu- 
ally it slowly vanishes in the last stage, which is dominated 
by radial fiow (IV, -Rmin = IQ-^-Ro). 



Fig. [3] displays the evolution of a symmetric, large ini- 
tial opening-angle shape, 6 = 75°, initialized with radial 
flow. Fig. [5]3 shows that the shape quickly evolves into a 
nearly-cylindrical segment connecting two cones. As the 
pinch-ofF proceeds and we zoom in close to the minimum, 
the shape becomes long and slender. The effective half 
opening-angle of the innermost profile in each sub-figure 
is (a) 6lcff = 75.0°, (b) O^s = 32.3°, (c) e^s = 13.3°, and 
(d) 6cff = 7.6° ^. This is a very fast process compared to 
the final stage, which roughly begins where (d) leaves off, 
where the shape very slowly approaches a cylinder. Fig. 
Ulshows why this takes place. Immediately after starting, 
the water develops a pressure peak just inside the neck. 
In order to satisfy the boundary condition P — Pair on 
the surface, the pressure peak takes on an aspect ratio 
similar to the neck aspect ratio. The pressure gradi- 
ents generate vertical velocity on the surface compara- 
ble to -Rmin- Although is zero at the neck minimum 
due to the symmetry, it becomes large just a short dis- 
tance along the surface. Fig. [5] shows the surface velocity 
fields. The inset shows that the vertical velocity diverges 
approximately in pace with the radial velocity. There- 
fore, we consider the rescaled vertical fiow, Uz {z, t) / RmXm 
in the neighborhood of the minimum (Fig. [SJ versus z 
rescaled by Lz). The rescaling shows that Uz always re- 
mains a small fraction of the radial collapse, reaching a 



maximum in the transient period (II) and then broaden- 
ing in spatial extent. 

We show the continuous evolution of the velocity distri- 
butions for 9 = 75° by plotting the maximum of Uz/Rmin 
versus the aspect ratio (Fig. [5^, uppermost curve). Ro- 
man numerals indicate corresponding stages between Fig. 
[5] and Fig. [6^. Just after t — 0, a stretching flow builds 
in strength before the shape has time to react (I). After 
reaching a peak vertical-to-radial ratio (II), the neck con- 
tinues to stretch (III) . The exterior flow decays slowly to 
radial implosion. Starting with different initial opening 
angles yields the same qualitative behavior. Smaller ini- 
tial opening angles generate smaller maximum Uz/ Rmin 
in the transient regime (I-II). Eventually all initial con- 
ditions enter the slender-body regime (IV) where we find 
the slow evolution towards a perfect cylinder described 
in previous works 




cff defined by 2R^ 



lin 



Figure 6: (Color online) Normalized maximum vertical ve- 
locity, max |uz I /-Rmin, versus shape aspect ratio, for sym- 
metric shapes with initially radial flow, labeled by 6 in de- 
grees (from top to bottom). The flow invariably becomes a 
stretching flow, increasing the aspect ratio of the break-up 
region, L^/Rmin- The normalized maximum velocity initially 
increases rapidly (I), peaks (II) and flnally decays onto one 
curve (III, IV). 

Discussion-While our analysis has focused on axisym- 
metric dynamics, the conclusions should remain relevant 
when the neck shape has a slight azimuthal asymmetry 
since prior linear stability studies have shown that the 
asymmetry does not grow rapidly but instead persists 
as vibrations of fixed amplitude. In the opposite regime 
when the neck shape is strongly distorted from axisym- 
metry, experiments show the bubble neck rips itself apart 
in a series of violent, apparently discrete jerks Q. The 
force balance controlling this break-up mode remains an 
open question. 

Conclusions-We showed here that strong, transient 
vertical flows are the dominant effect in the transition 
of the neck shape from large to small cone angles. This 
transition happens as a distinct, transient phase, as op- 
posed to the very slow, continuous transition predicted 
by previous studies assuming weak vertical coupling. The 
transition is effectively complete once the minimum ra- 
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dius has decreased by a factor of 100, an easily observable 
window. Even if the initial cone angles are very large, the 
system just induces a stronger vertical flow, still guaran- 
teeing a rapid transition. Moreover, if the vertical veloc- 
ity is initially compressing, it quickly reorients to a strong 
stretching flow. All these things taken together indicate 
that the transition to small cone angles will occur before 
any cut-off length scale. After the magnitude of the verti- 
cal flow peaks relative to the radial flow, it soon becomes 
small in comparison. This guarantees that the system 
will reach the slender-body phase characterize d by a log- 
arithmically slow transition to cylindrical Ill l 4l3l. along 
with memory of azimuthal vibrations [tI. Il6l]l7|. both 
studied extensively. In sum, the radial-flow dominated 
implosion singularity controls the final stage of the dy- 



namics even in situations with strong vertical coupling. 
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